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Abstract. We use SDSS spectra and optical to far-infrared photometty for a sam- 
ple of 3 1 FeLoBAL QSOs to study the relationship between the AGN-driven outflows, 
and obscured star formation in their host galaxies. We find that FeLoBAL QSOs in- 
variably have IR luminosities exceeding 1O 12 L . The AGN supplies 75% of the total 
IR emission, on average, but with a range from 20% to 100%. We find a clear anti- 
correlation between the strength of the AGN-driven outflows and the contribution from 
star formation to the total IR luminosity, with a much higher chance of seeing a star- 
burst contribution in excess of 25% in systems with weak outflows than in systems 
with strong outflows. Moreover, we find no evidence that this effect is driven by the 
IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN act 
to curtail obscured star formation in the host galaxies of reddened QSOs to less than 
~ 25% of the total IR luminosity. This is the most direct evidence yet obtained for 
'quasar mode' AGN feedback. 



1. Introduction 



Over the last decade, several problems have arisen in modelling the cosmological evo- 
lution of galaxies. These include (1) the difficulties that models faced in explaining 
the observed galaxy luminosity function at low and high redshifts simultaneously (e.g. 
Benson et al.1 120031) . (2) the prediction t hat rich galaxy clust ers should harbour cool- 
ing flows when few are observed (e.g. IPeterson et al.l 120031) . and (3) the difficulties 
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that models face in reproducing the number of IR-luminous galaxies observed at high 
redshift (e.g. iBaugh et alj|2005h . 

One of the most promising solutions to these issues is AGN feedback. AGN feed- 
back is the exertion of influence of an SMBH on >kpc scales to curtail star formation 
in the host galaxy, and/or accretion onto the SMBH itself. Galaxy evolution models 
usually assume this feedback to occur in one or both of two simplified modes; 'quasar' 
mode and 'radio' mode. Quasar mode feedback occurs via radiation from an accretion 
disk, while radio mode feedback occurs via a relativistic jet that transfers momentum 
to the ISM. Both these feedback paradi gms have led to improvements in the abil ity of 
models to reproduce observations (e.g. Bower et"al]|2008l : ISomerviile et al.ll2008h . but 
observational evidence for either paradigm remains sparse. 

Our group has been looking for evidence of quasar mode feedback by examining 
systems in which s uch feedback may be ongoing. To do so, we have been studying the 
'FeLoBAL' QSOs dHazard et all 1 19871 : lHall et al.ll2002l) . We selected this population 
for two reasons. First, their UV absorption troughs are unambiguous signatures of ra- 
diatively driven outflows powered by an A GN. Second, FeLoBAL Q SOs are invariably 
reddened and have high IR luminosities dFarrah et al.ll2007l l2010h . We here present 
results for 31 FeLoBAL QSOs, in which we compare the strength of their outflows as 
estimated from their UV spectral properties, to the luminosity of obscured star forma- 
tion in their host galaxies as measured from optical through far-IR photometry. We 
define "IR luminosity" as the luminosity integrated over 1-1000/mi in the rest-frame. 



2. Methods 



We assembled our sample from the SDSS. We started with the six SDSS objects in 
Farra h et al.l2007l (hereafter F07) and added a further 25 FeLoBAL QSOs from lTrump et al. 



(2006) over the redshift range 0.8 < z < 1.8 (chosen so that the SDSS spectra always 
contain the Mg II BAL). There was no selection on IR luminosity. We observed the 
sample with MIPS ( Ri eke et all 120041) onboard Spitzer using standard parameters. We 
reduced the MIPS data using the MOPE X software. We then added archival photom - 
etry from t he SDSS, 2MASS o r UKIDSS (ISkrutskie et alJl2006t Eawrence et al.ll2007l) 
and WISE dWright et al.ll2010r) to give photometric coverage from 0.9//m to 160/im 

We measure total, starburst and AGN IR luminosities by fitting the optical through 
MIPS photometry for eac h object with radiat ive transfer models for AGN and starbursts, 
following the methods i n Farrah et al.ll2003l and F07. For the AGN models, we follow 
Efstathiou & Rowan-Robinson ( 1995b and u se a tapered disk dust distribution. For 



the starburst models we use the templates o f lEfstathiou et all ( 2000!) with the updated 
dust model of lEfstathiou & Siebenmorgerj ((2009). The models do not include broad 
absorption features in the rest-frame UV that arise from radiatively driven outflows, so 
we do not include the SDSS spectra, or photometry shortward of ~0.35/mi in the rest- 
frame, in the fits. In 27/31 cases the models provide a good fit to the data (jC r d < 2). 
The remaining four objects have 2 < x 2 . < 3, usually due to a poor fit in the near-IR. 
This could be due to the lack of a host galaxy contribution in our models, so we defer 
exploration of this until more comprehensive near-IR data are available. 
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Figure 1. Absorption strength vs (left) starburst (stars) and AGN (squares) IR 
luminosities, and (right) starburst contribution to the total IR luminosity. 



3. Results & Discussion 

We find that FeLoBAL QSOs are luminous in the IR. All 31 of our sample have to- 
tal IR luminosities (Lrat) in excess of 10 12 L o , with nine objects exceeding 10 13 L o . 
These luminosities ar e comparab le to those of the wider population of BAL QSOs 
dGallagher et al.ll2007h . red QSOs dGeorgakakis et al.ll2009h . and ULIRGs dFarrah et al. 
120031. 120091) . The dominant power source is usually an AGN. A pure AGN is either the 
most likely power source, or consistent within the 90% confidence interval, for 35% of 
the sample. A starburst component is required for the remaining objects, but in only 
twelve objects is the starburst more luminous than 10 12 L o , and in only three objects 
is the starburst more luminous than the AGN. The mean starburst contribution to the 
total IR luminosity (fsg) is ~ 24%, compara ble to that of local U LIRGs with 'warm' 
IR colours, but lower than those of PG QSOs dVeilleux et al.ll2009l) . 

We now examine whether or not there is a relationship between the AGN-driven 
outflows and the obscured star formation. To quantify the strength of the outflows we 
use the same, single specie s (Mg II/12799A) across the whole sample, and adopt the Bio 
parameter of iGibson et al. (200 9J). We meas ured the Mgll absorption strengths using 
the methods described in lUrrutia et al.1 (120091) . 

We first examine if absorption strength depends on AGN or starburst IR lumi- 
nosity (Fig [TJ left). There is no correlation between absorption strength and L5 # (a 
Spearman rank test gives p = -0.10, P = 0.58) and at best a weak correlation between 
absorption strength and Lagn (p - 0.39, P = 0.04). Conversely, we see a weak but 
clear anticorrelation between absorption strength and f$B (p = -0.49, P = 0.005, FigQ] 
right). Moreover, all the systems with Bio 3500km s" 1 have fsB < 25%, while the 
systems with Bio ^ 3500km s" 1 have a wide dispersion in starburst contributions, from 
0% to ~ 80%. A two-sided Kolmogorov-Smirnov test for the objects above and below 
Bio - 3500km s -1 yields a difference in distributions at 99.84% confidence, though the 
number of objects in the two subsamples is too low for this test to be robust. 

We examine this result in more detail by constructing probability distribution func- 
tions (using all the fit solutions) for f$ b for two subsamples divided at Bio = 3500km 
s _1 (Figure |2j left). We see a clear difference. The low absorption strength subsample 
shows a much higher chance of a higher starburst contribution than the high absorp- 
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tion strength subsample. We quantify this by extracting the probabilities of obtaining 
fsB > 25%. For the whole sample we find P(f SB > 25%)= 50.3+^%, for the low 
absorption strength sample we find P(/sb > 25%)= 67.3^%, while for the high ab- 
sorption strength sample we find P(fss > 25%)= 17.8+H%. 

This anticorrelation between absorption strength and f$ b is straightforwardly in- 
terpreted as the outflow from the AGN curtailing star formation in the host galaxy. In 
this scenario, the systems with fss > 25% are those in which an outflow has yet to 
curtail star formation, so the outflows always have Bio ^ 3500km s -1 . Conversely, the 
systems with /sb < 25% are those in which an outflow has curtailed star formation, 
and those in which such an outflow has subsequently waned, making the dispersion in 
absorption strengths wide. 

There are however four other ways that this anticorrelation could arise. The first is 
that the starburst 'suppresses' AGN outflows, so when the starburst wanes (via a cause 
unrelated to the AGN) an AGN driven outflow can appear, possibly because the ISM 
density has been reduced by the starburst. We cannot test this scenario, but it would 
likely require a serendipitous conjunction of ISM and SMBH parameters, so we do 
not consider it further. The second is if a high starburst contribution meant that the 
Mg II troughs were observed to be weaker than they really are. Again, we cannot test 
this, but the rest-frame UV continua of starbursts in ULIRGs are at least an order of 
magnitude too wea k to provide this effect, and can show absorption in the same species 
(lFarrahetalJl2005h . The third is if BAL QSOs with strong starbursts preferentially 
drop out of the SDSS QSO s election compared to BAL QSOs with weak starbursts, 
and were thus not included in lTrump et al.1 (12006). This possibility is also not testable 
by us, but the SDSS is now turning up FeLoBAL QSOs in large numbers, and the SDSS 
QSO followup colour selections are fairly relaxed, so we do not consider this possibility 
likely either. 

The fourth possibility is that stronger outflows reflect an increase in the IR emis- 
sion from the AGN, but have no effect on the starburst; if this is the case we would 
see a decline in the contribution from star formation to the total IR luminosity with 
increasing absorption strength, but with no direct relationship behind the decline. This 
is a possibility we can test, as follows. If outflow strength is just a proxy for AGN 
luminosity, then we should see a bigger difference between the starburst contribution 
PDFs for subsamples divided by AGN luminosity than between subsamples divided by 
absorption strength. In the right panel of Figure |2] we show starburst contribution PDFs 
for two subsamples divided by AGN luminosity at Lagn = 1O 125 L . The difference 
between the PDFs divided by AGN luminosity is weaker than the difference between 
the PDFs divided by absorption strength. Furthermore, we are more likely (albeit only 
at ~2<t) to obtain a smaller starburst contribution by selecting high absorption strength 
systems than we are by selecting high AGN luminosity systems; for P(/sb > 25%): 
the BI > 3500km s" 1 subsample is 17.8+^% while the Lagn > 10 125 L o subsam- 
ple is 38.7^ 5 Q %). In other words, we are more successful in finding systems with 
a large starburst contribution to the total IR emission by selecting on weak outflows 
than we are by selecting on low AGN luminosity. Also, we are more likely to find star 
formation with a lower absolute luminosity in the lower luminosity AGN subsample 
(e.g. for P(L sb < 10 12 L o ): the Lagn < 10 125 L G subsample is 69.4^- 2 % while the 
Lagn > 10 12 5 L Q subsample is 43.2+^%) 

Overall therefore, we find that radiatively driven outflows from an AGN with ab- 
sorption strengths > 3500 km s -1 act to curtail star formation in their host galaxies. 
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Figure 2. Probability Distribution Functions for the starburst contribution to the 
total IR luminosity. In both panels, the solid line is the PDF for the whole sample. 
The other lines are for subsamples divided according to two criteria. Left panel: by 
absorption strength. Right panel: by AGN IR luminosity. 



We also find that this effect is (at least largely) relative; such outflows reduce the con- 
tribution from star formation to the total IR luminosity to less than ~ 25%. We also 
propose that the infrared luminosity of the AGN is not a good proxy for the degree of 
AGN feedback that is taking place. 
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